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Abstract 

The Cantor set complementary to the Devil's Staircase associated 
with the Circle Map has a fractal dimension d <~ 0.87, universal for a 
wide range of maps, such results being of a numerical character. In 
this paper we deduce a formula for such dimensional value, the corre- 
sponding theoretical reasoning permits conjecturing on the nature of 
its universality. The Devil's Staircase associated with the Circle Map 
is a function that transforms horizontal unit interval / onto /, and is 
endowed with the Farey-Brocot (F — B) structure in the vertical axis 
via the rational heights of stability intervals. The underlying Cantor- 
dust fractal set ft in the horizontal axis, ft C I, with fractal dimension 
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d(Cl) = d <~ 0.87 has a natural covering with segments that also fol- 
low the F — B hierarchy: the staircase associates vertical I (of unit 
dimension) with horizontal Cl C / (of dimension ~ 0.87), i.e. it selects 
a certain subset of /, both sets F — B structured, Q with smaller 
dimension than that of /. Hence, the structure of the staircase mir- 
rors the F — B hierarchy. In this paper we consider the subset Qf-b 
of I that concentrates the measure induced by the F — B partition 
and calculate its Hausdorff dimension, i.e. the entropic or informa- 
tion dimension of the F — B measure, and show that it coincides with 
d(fl) ~ 0.87. Hence, this dimensional value stems from the F—B struc- 
ture, and we draw conclusions and conjectures from this fact. Finally, 
we calculate the statistical "Euclidean" dimension (based on the ordi- 
nary Lebesgue measure) of the F — B partition, and we show that it is 
the same as d^lp-B), which permits conjecturing on the universality 
of the dimensional value d ~ 0.87. 

1 Introduction 

The Cantor set complementary to the Devil's Staircase associated with the 
Circle Map has a fractal dimension d ~ 0.87 [Jensen et al., 1984], univer- 
sal for a wide range of maps [Bak, 1986], such results being of a numerical 
character. In this paper we deduce a formula for such dimensional value, 
the corresponding theoretical reasoning permits conjecturing on the nature 
of its universality. 
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A Cantor or Devil's staircase is an increasing function y = g(x) from 
/ = [0, 1] onto /, with zero derivative almost everywhere, constant in the 
so-called intervals of resonance or stability Axj,, k 6 IN, which are infinite in 
number. Such staircases are frequently observed in empirical physics [Bak, 
1986], and their universal properties are of great interest. The complement 
in / of Uk<=wAxk is a totally discontinuous Cantor-dust set f2 naturally asso- 
ciated with the staircase, which reflects the features of the particular physical 
problem under study. The sine circle map n+ \ = 9 n + w + ^ sin(27r# n ) is 
a simple model describing [Bak, 1986] systems with two competing frequen- 
cies, e.g. the forced pendulum, with 9 the angle formed by the vertical and 
the pendulum, n the discretized time variable, and w the frequency of the 
system in the absence of the non-linear term given by the sine function. Let 
W = lirrin^oo ^ be the winding number of the system. The graph of the 
function W = g(w) is a well known Cantor staircase; with Aw we denote 
an interval of stability as well as the corresponding stair step. 

Let Awi and Au)2 be two such resonance intervals such that all intervals 
between these two have smaller length. Let g(x) = | if x G Aiui, g(x) = y 
if x G Ati>2, all stair steps have rational height. If x is in the largest interval 
in the gap between Aw\ and A^2, then g(x) = fqEp-; i.e. the height of stair 
steps follows what, by definition, is the Farey-Brocot (F — B) interpolation 
law. This is so for many staircases empirically found in physics and other 
sciences. Starting from | = ° = 0, and y = \ = 1, F — B interpolates 
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= \ between and 1, partitioning / in two intervals, in turn partitioned 
in two intervals each, yielding a partition of / in 4 intervals in the second 
order interpolation. . . and in 2 N intervals in the N th order of the F — B 
interpolation. The induced F — B measure in level N of interpolation gives 
the same probability measure — by definition — i.e. to each of these 2 N 
intervals. Let Cm = I — Up Awe be such that wv is the stair step of height 
|, and | is every rational in the N th level of F — B interpolation. Then Cm 
is a covering of £1 by 2^ intervals Ij N ^ in the horizontal w axis, such that 
g(Ij N ^), j : 1 —> 2 N are the 2^ intervals of the F — B partition of / in the 
vertical W axis. If we plot [Piacquadio, 2004] lengths of the 2^ intervals 
Ij N ^ against length of the corresponding g(Ij N ^) we obtain a straight line 
that passes through the origin with slope cjy growing as N grows. So 
also follows the hierarchy of the Farey tree via its covering: the staircase 
relates an F — B structured unit segment / (dimension 1) with an F — B 
structured subset (dimension ~ 0.870 ) of /, i.e. the staircase selects a 
subset of / of smaller dimension, the F — B partition being at the core of 
the very structure of the staircase. Hence, it seems natural to relate the 
F — B partition to the ~ 0.870 value: using the tools of multifractality, we 
calculate the multifractal spectrum (a,f(a)) of the F — B measure on /, 
and identify which subfractal Qf-b C / has a dimension ~ 0.870. We find 
that Qf~b is the set that concentrates the F — B measure, it corresponds 
to the value a for which f(a) = a and f'(a) = 1, i.e. its dimension ~ 0.870 
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is the entropic or information dimension of the F — B measure. 

A feature in the importance of this dimensional value is its universality, 
which has been checked (again [Bak, 1986]) by studying a broad class of 
circle maps with more complicated non-linear terms than the simple sine 
map. Although details may differ from those of the sine map W = g(w) 
— steps narrower, sometimes larger — still the dimension of the underlying 
remains ~ 0.870. 

We proceed as follows: the 2 N intervals in the N th F — B partition have 
the same F — B measure, ^r, but very different lengths, i.e. very different 
Euclidean measure. We start (Sees. 3 and 4) with the thermodynamical 
algorithm (a, f(a)) for the hereinafter called Euclidean case, and by this we 
mean: all segments considered have equal Euclidean length at any N th parti- 
tion. We proceed from there in slow steps in such a way that the results can 
be extended to the F — B measure (Sec. 5) in a manner that — we trust — 
will be seen as "natural" . Thus, a first connection between the two measures 
will be established: working always in /, we express a and f(a) in terms 
of contractors (probability contractions pj G (0,1) and/or length contrac- 
tions Cj £ (0, 1)) and their key frequencies linked to each other through the 
thermodynamical algorithm; a finite number of contractors for the so-called 
Euclidean case, extending the results to an infinite number of contractors 
in the F — B case. Next, we estimate (Sec. 5.5) the Hausdorff dimension of 
the subfractal Q.F-B for which f(a) = a and f'{a) = 1 for the F — B mea- 
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sure, and we obtain the entropic or information dimension of the measure to 
be -5 Aii lgg^0 ± ^' wi th A, = £ and c = which yields the 

value 0.87038 in the interval defined by 0.870 ±0.0004 (again Bak 1986), the 
universal constant associated with the dynamics of the Circle Map — which 
is why we conjecture that said dynamics inherits, via the F — B structured 
staircase, this universal constant, which is an inherent property of the F — B 
measure. 

Finally, by taking averages over the very different lengths of intervals in 
Nth F — B partitions, as N grows, we obtain a statistical F — B contractor 
(Sec. 7) with which we can calculate the dimensionally Euclidean (having 
only one contractor, all segments have equal Euclidean length at any Nth 
partition), statistically self-similiar, fractal version of the F — B partition. 
Such process, briefly described in Sec. 7 yields again the universal value 
~ 0.870 (Sec. 8), which is a second and much deeper connection between 
the two measures. 

The definition of the F — B measure on I as constant over intervals in 
Nth partitions is not an arbitrary one: there is a non-Euclidean geometry 
on the upper half-plane, the F — B partition is its inheritance on /. This ge- 
ometry (we briefly comment on it in Sec. 9) has an associated regular tiling 
which partitions the real line in F — B interpolations, and the location of 
the tiles approaching a real irrational number " i" describes — by naked-eye 
direct observation — its decomposition in continued fractions, which yields, 
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as we will see, the location of " V in the multifractal spectrum of the F — B 
measure. 

NOTE: Sec. 6 can be by-passed by the reader: there are old and new 
results on the (a, f(a)) spectrum of the F — B measure, and in Sec. 6 we 
"harmonize" the — only apparent — corresponding discrepancies. 

NOTE: Although the Math level in this paper does not go beyond find- 
ing extremes of a function of several variables, the reader un-interested in 
long and tedious and tiresome estimates, approximations, and calculations 
can proceed to Sec. 2: Generalities and Notations, then go to Sec. 7 and 
therefrom to Sec. 9: Geometrical Considerations, Conclusions and Conjec- 
tures. 

IMPORTANT NOTE: By necessity we work, from Sec. 3 to Sec. 8, 
sometimes with approximations " =" , sometimes with exact equalities " =" , 
so the corresponding calculations yield estimates, and make no claim to be 
rigorous proofs of formal theorems. 

2 Generalities and Notation 

With pj we will denote probabilities, A-,- will be frequencies, Cj contractors, £ 
will be a "normalizing sum" (a different one for each normalizing process); A 
and fi the coefficients of the Lagrange method of indeterminate coefficients 
for finding extremes of functions. With "i" we will denote an irrational 
number in the unit interval, i = — ■ 1 i := [a±, a,2, ■ ■ ■ a n , • • •] is its contin- 
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ued fraction expansion; a,j the so-called partial quotient coefficients will be 
natural numbers; the rational number \a\ . . . a n ] := — is the nth rational 
approximant to i; q n the so called nth cumulant. 

The Farey-Brocot tree interpolates rational between rationals | 
and y, starting from ° and j, the extremes of the unit interval. The 
first interpolation has, therefore, two segments [j, |] and [|, y], . . . the Nth 
step of interpolation partitions the unit segment in 2^ segments. Approx- 
imant ^p- = \a\ . . . a n ] appears in the Nth step of the partition process, 

N = oi + + a n ; i G ^+1] Vn G iV (provided that 2™ < £n+i) and 

the length of this segment is q ^ — . For a certain probability measure on 
the unit segment, let us consider a partition, lj the length of its segments, pj 
their probability measure. Let q and r, real numbers, be connected through 
Y2j jrr = 1, r = r(g). Then, the so called thermodynamical formalism yields 
the multifractal spectrum (a,f(a)) of the probability measure in terms of 
lagrangian coordinates: a = r'(q)\ f(a) = aq — r; q = f'(a); f"(a) < 0, and 
q = 1 iff a = /(a). Let us recall that f(a) is, theoretically, the Hausdorff 
dimension of the subfractal O a which contains all elements with the same 
a-concentration; and that the a-concentration of a segment is the log-log 
version of the density: ^||, I the length of the segment and p its probability 
measure. Point a-concentration is defined in the same way as point den- 
sity: from a-concentration of segments containing the point and a limiting 
process. 
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3 The Euclidean Case: Equal Lengths 
3.1 Two initial probabilities 

First we consider the unit segment as partitioned in two segments of equal 
Euclidean length when N = 1,. . . and in 2^ Euclidean equal parts in step N. 
The first two segments have probability measures p\ and p2 , p\ + P2 = 1 , so 
Pj , j = 1,2 are contractions. The four segments of equal length in step N = 2 
have probability measures pi 2 , P1P2, P2P1 and P2 2 respectively . . . and so on. 
In step N, we have 2 N segments, and their generic probability is pi r p2 N ~ r , 
r an integer, < r < N. The number of segments with this probability is 
{ N r ) ■ We will redo a calculation of the multifractal spectrum (a, f{a)) of this 
measure in terms of a key frequency A = as internal coordinate. In step N 
we have, for an "r" segment, a concentration a = lo ZP rob - = lQ g(pi Pg — ) _ 

' to ' loglength iff(pr) 

Ho gPl _ + ^)lo gP2 = _=A_ {j _ logpi + (1 _ r_ )logp2) = ^(Alogp! + (1 - 

A)logp2)- We can write then a(A) for this magnitude. The number of 
such elements is (^) = , t ^ 1 -t, = r ,^\ N -r \ = H r — 1 — ~} N = 

\ r ) r\\N—r)\ r r (N—r)" r e e L /" JL~i7v n ?Ly~iv J 

\{ \\(x-\y->> ) N • Therefore, in that step, and for such a, we have f(a) = 

e A i g 2 y A = Ep{ A1 °g A + (! " A ) lo g(l ~ A )i if N tends to °°- 0ur 
system now reads a = ^(Alogpi + (1 - A)logp 2 ) = 15^ (A log p + (1 - 

A)log(l-p)); pi :=p; f(a) = To ^ 2 (AlogA+ (1 - A)log(l - A)), from which 



A 

-A 



it is obvious that a = /(a) iff A = p. Also /'(a) = £ = |^ = 
and again, /'(a) = liffA = piffa = f(a). As variable A : — > 1 varies, the 
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(a, /(a)) graph is drawn. We call A an "internal" coordinate, because in a 
certain step N, A = tells us the value of r, the number of times in which p 
appears, that is the proportion ^ of all segments with measure p\P2 N ~ r = 
p r (l—p) N ~ r . Henceforth, we will be interested in the subfractal for which 
a = f(a) and f'(a) = 1, in all cases and in all measures considered. That 
we have, for this case, f'(a) = 1 iff a = f(a) validates the thermodynamical 
formalism (see Sec. 2), which is not at all proved to yield the multifractal 
spectrum of an arbitrary measure, but which holds true for the Euclidean 
measure. In this case Ylj frr = 1 for /'(a) = g = 1 yields Ylj = 1 
which, together with Y^Pj = 1) an d the fact that partition {lj} is arbitrary, 
imply r = 0, which is another feature of the thermodynamical algorithm: 
a = /(a) iff q = /'(a) = 1 iff r(q) = 0. 

Next, let us arrive at the expressions for a, f(a), and /'(a) through 
the lagrangian coordinates in the thermodynamical algorithm, and compare 
said expessions with those above, with A : — ► 1 as the internal coordinate. 
From Ylj = 1) an d h = jn V?', in step N, we have 

„ logEjPj* 1 . vv 17. i»-n«W 

T= ^- = -m^2 l °^^-p) y[ r ) = 

— 1 N fN\ — 1 

■ log^WKi -p) q ] N - r ( = m^ l °s(p q + (i-p) q ) N = 

, — n \ / o 



iVlog2 



1 log(p' + (l-p)«). 



log 2 
So 

a = r'(q) = — — 7^ r-ia 13 logo + (1 — p) q log(l — p)\ = 

w log2p<? + (1 -p)i l V 7 V n 
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1 r P q i , 



log 2 V + (1 -p)i + 
a comparison with a = [^{Alogp + (1 — A) log(l — p)} yields 

A= f- (1) 

pi + {l-p)i y ' 

Next, following the algorithm, we have f(a) = aq — r = 

-1 p q (1-p) 9 -1 



log2 l pi + (l-p)i & ^ £ m log2 
-1 (1 - v) q V q (I - v) q 



log2 L £ °^ £ m v ° /L £ £ 
-{ — log — + log } 



log2 £ £ £ ° £ 

and again, if we compare with the value [^{AlogA + (1 — A)log(l — A)} 
we obtain A = ^. Now, /'(a) = |° g 1 ^ P X , with this value of A, and with 
= ' srnce £ cancels, becomes 



indeed. 

This will be the procedure for the next sections: to express a and f(a) 
in terms of contractors and key frequencies, to find an expression for these 
frequencies in terms of the thermodynamical parameters q and r, and to 
find the frequencies for which a = /(a), /'(a) = 1, and r = 0. 



3.2 A finite number of initial probabilities 

Consider, next, the unit segment as partitioned in no segments of equal 
Euclidean length when N = 1, . . . and in riQ N equal parts in step N. The first 
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no segments for N = 1 have probability measures p\ , . . . , p no ; pi + . . . +p no = 
1; Pj contractors. This case is quite different from that in which uq = 2: 
the frequencies of p\ and p2 = 1 — p\ were Ai and A2 = 1 — Ai, so there 
was a coordinate A = — > 1, a "natural" or "internal" coordinate in charge 
of producing the spectrum. We cannot have that convenience here, for the 
frequencies of the pj, the Xj, will be Ai . . . A no , 2^ Aj = 1, so we have many- 
independent coordinates. 

Let N be the step, r\ . . . r no a particular choice of integers, < rj < N; 

Tj = N, we consider segments of length -\ with measure p'^ . . . p^ ° ■ We 
proceed as in the previous section: Xj = ^, Xj = 1, const • — ^7 — ^ the 
number of such segments with the {rj} or the {Xj} particular choice. The a- 
concentration and the f(a) corresponding to such a set {Xj} are: a(X\ . . .) = 

A -°iognt ' •••)) = Al l°iog»t '" ' P roceedin g as in the previous section. 

But the difference with last section arises now: for a fixed value of a we are 
interested in all choices of {A^} which fulfill a(X\ . . . Xj . . .) = a. And the 
dimension f(a) of this subfractal will be the maximum value of f(a(Xi . . .)) 
which fulfills a(Ai . . .) = a and ^ Xj = 1. Therefore, we have to extremize 
-I^(Ai 1 ogAi + ...)-A I ^(A 1 logp 1 + ...-a) + / u(Ai + ...-l), with Xj 
as variable. The corresponding calculations are shown in the App. to Sec. 
3.2; the result: 



A, = f Vj (3) 
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with 

A = A(a), or a = a(A) (4) 

which yields /'(a) = A. 

So, the lagrangian indeterminate coefficient A fulfills four roles: (1) it 
is the lagrangian coefficient linking /(a) with a; (2) it is the exponent, in 



24 



Eq. 23, of pj, which, normalized, determines A^; (3) it gives a from Eq. 
and (4) it is f'{a) = q, as we have just seen. Notice the similitude between 
these results and Eqs. m ^1 From the values of a and /(a) obtained in the 



Appendix to Section 3.2 (see Eq. 25), we can see that a = f(a) iff \j = pj, 



which happens iff A = 1, i.e. if f'(a) = 1. All of which implies r = 0. 



4 The Euclidean Case: Equal Probabilities 

In this case the lengths of all Uq segments in the partition of the unit segment 
corresponding to the Nth stage or step of the construction of the multifractal 
are given by contractors c\ . . . c no , a natural extension of the case in Sec. 3. 
All of the n,Q probabilities are equal. With Xj as before, the generic length 
of such a segment is (c^ 1 . . . Cno°) N ■ Proceeding as in Sec. 3. we have to 
extremize the function ^gg^ - A( a^hL ~ a ) + M^i + . . . - 1) with 
Xj as variable. The corresponding calculations are shown in Appendix 1 to 
Sec. 4; the result: 
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which bears a resemblance to Eq. |3j here 3 is an exponent of the contractor 
Cj, and 

f («) = i (6 

log n 

We want to interpret Eq. [6j since it does not look anything like. 

Let q = f'(a),T(q),a,f(a), be the thermodynamical magnitudes in- 
volved in the process described in Sec. 3: equal lengths and different proba- 
bilities. In Sec. 4 we are reversing the process, exchanging the role of lengths 
and probabilities: equal probabilities and different lengths, which has been 
termed "the inverse process" . Notice that ^ • pj = 1 and ^ ■ Cj = 1 make 
this inversion totally plausible. Let /, a, q, f , be the new thermodynamical 
parameters. 

In Appendix 2 to Section 4 we deduce the relationships between "old" 
a i fiQ = f'( a )i T an d the new /, a, q, f: 

a = -, (7) 

- r(q) = q Vg (8) 
-f(-r(q)) = q Vg or - r(-r(g)) = q Vg (9) 

which yield 

q= — = - =-r(a) (10 

da log no 

and 3 = — f(A), from which q = — r(— f(A)) = A or q = A in Sec. 4, as 
q = A in Sec. 3. 

The expression 3 = — f (A) above becomes 3 = —f(q), which is q, by Eq. 
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29 — that is, the "old" q. Hence, cases in Sees. 3 and 4 have an analogy: the 
lagrangian coefficient A is /' in both spectra, and a difference: the exponent 
of the contractor giving the critical A is, in the first case, A = q and, in the 
second case, the q of the inverse problem; i.e. the A of the inverse problem. 

Notice that condition a = f{a) (in Sec. 3/Sec. 4 notation we should 
write a = f(a)) is fulfilled for A, = — , for then Ai log Ai + . . . = — log — + 

■i \ 1 1 j no' no ° no 

. . . = no^(— logno) = —log no, the numerator of a. But Aj = which 



from Eq. 26 means S = 0, i.e. — r(A) = 0, which means A = f'{a) = 1. 

Again f{a) = a, f'{a) = 1, r = are simultaneous conditions in order 
to characterize the subfractal which concentrates the measure, /(a) being 
the entropic or information dimension. 



5 The F — B Case 

5.1 Equal probabilities and different lengths 

The treatment of the thermodynamical multifractal spectra in the Euclidean 
case, expressing key parameters in terms of contractors and their frequencies 
in Sees. 3 and 4, permits — we trust — extending such results and reasonings 
to the case of the F — B measure on the unit interval given by the Farey 
Brocot F-B partition tree. As in Sec. 4 we deal with equal probabilities and 
different lengths. 

The Nth step or stage of the F-B interpolation gives a partition of the 
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unit segment in 2^ smaller segments of equal t^v probability. Let a\ . . . a n 
be positive integers such that a\ + . . . + a n = N . There is a segment in 
that step of length — - — , where \a\ . . . a n ] = ^ (see Sec. 2). This segment 

Qn—lQn L 9n 

contains all irrationals of the form i = [ax, . . . ,a n , etc. ], where "etc." is any 
sequence of natural numbers a n+ \, . . . , a,j, . . .. 

We want to interpret nested segments of length — - — , — - — , . . . in 
terms of contractors. 

5.2 The F — B contractors 

First we observe that lengths of nested intervals diminish like \: since 



Qn+l 



a n +iq n + q n -i, the q n grow with n, hence -J— < qn g n+1 < ^i- 



Therefore, we will estimate lengths — — - — , q n = q n (ai, . . . , a n ) by 4- in step 
N = a\ + . . . + a n . Now, 

a n +iq n < q n +i = a n+ \q n + q n -i < a n+ iq n + q n = (a n+ i + l)q n . (11) 



So, the contractor that shrinks length -4- into the smaller one -4— is a 
number somewhere between -4— and -, 1 . Now, a moment of reflection 



observing Eq. 11 shows that q n +i is much nearer (a n +i + l)q n than a n +iq n : 



o-n+iqn is far smaller than q n +i because an exponential (q n -i) is missing, 



whereas by replacing q n -\ by q n (in the RHS of Eq. 11) we just replace 
one exponential by another which could be connected to the first one by a 
reasonable coefficient. 

Therefore, a certain contraction c(a n+ \ + 1) of a n+ \ + 1 will yield q n +i 
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from o„, via Eq. 11 hence -J— from \. 

Let us consider q n = g„(ai . . . a n ). Magnitude q n is obtained from q n -\ 
via c(a n + 1), we still do not know the value of c. Iterating this process we 
have q n (cL\ ■ ■ ■ a n ) given by c{a\ + l)c(a2 + 1) . . . c(a n + 1). The n integers 
a\ . . . a n vary — say — between 1 &; k G IN. Let r,- be the number of times 
for which the a's are equal to j : 1 — * k, Y2j r j = n - Then q n is given by 
[c(l + l)] ri . . . [c(k+l)] r ' k = c ri+ -+^2 ri . . . (k + l) rk = c n [2 Al . . . (k + l) Xk ] n = 
[c2 Al . . . (k + l) Afe ]™. Here A; is, simply, the largest of the integers a\ . . . a n , 
Xj = and Ai + . . . = 1. 

We rewind: segment of length is obtained from that of length \ 

1n+l In 

through a contraction — i a n+i an integer. The estimate 

q n ^[c2 x K..(k + l) x T (12) 

above, c an appropriate constant, is a simplified version of the Besicovitch 
formula [Good, 1941], which we have already used elsewhere [Piacquadio, 
2004]. We are in the F-B step N = a\ + . . . + a n , in a segment of length 
estimated by \ = 7— r — ,} , k simply the largest value of the a,-, 

j : 1 — > n. The probability (F—B) measure of such segment is ^ 



2^ _ 2 E "=l^' 

We need, now, to estimate constant c € (0, 1). This we do in the Ap- 
pendix to Sec. 5.2, by estimating the Hausdorff dimension dn of E^ := {i = 
[a\ . . . a,j . . .]/a,j < k Vj}, k £ IN, Aj the frequency in which the a's are equal 
to j. We use a result of Jarnik [1928; 1929] who proved 1 - d H {E k ) = 0{\), 
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and obtain the Xj responsible for the dimension: 

A, = V±_L (13) 

with 

a = -2d H (E k ) (14) 

a = —2 as k grows, which implies c = ^/c^; c n = ^ — 1. 

From our "We rewind" note above, we have the generic value of the 
F - B contractors: [^fc^+ip and , since a n+i is any integer j, 
is the generic contractor, j G IN. The main difference with the Euclidean 
case is that we have an infinity of contractors now. 

5.3 A first estimate of the (a,f(a)) spectrum for the F — B 
measure 

From the preceding section the probability of segment with length l/[c2 Al . . . [k+ 
l) Afc ] 2n is 1/2 N = 1/2^" a i\ k the largest of the a's. The a-concentration of 
this segment is, then, 

a = -log2 = log 2 (Yli a i)/ n . = ]og2 m 

-2n logc + Ailog2 + . . . 2 log c + Ai log 2 + . . . ' 2 log c + Ai log 2 + . . . 

(15) 

Integer n "disappears" in the average value m of the a's, whereas k will 
become quite relevant. 

Let us consider, in E^, the set S m of elements with average of the a's no 
larger than m — technically, it should be lim sup n (^" a.j)/n < m, but the 
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essential idea is to control the average of the a's. A choice of A's: {Ai . . . A&} 
will label different subsets of S m . As we saw above, the subset of largest 
dimension corresponds to the label A 7 - = — , and we have to add the 
extra condition on the size of the average of the a's. This particular choice 
of A's is both responsible for the dimension of S m and, therefore, for the 



value of a associated with it, which, from Eq. 15 becomes a = t£,K the 



denominator in Eq. 15 for these particular A's. 



Now, let Vj be, as before, the number of a's equal to j, Xj = ^r, then 
= i Al + 2A 2 + . . . + k\ k = ± YX (jflF = ^ < m ...so k cannot 
be larger than k m = e c ^ m = e C7raK := e Ba , tl B" a constant. Applying 
the already quoted result by Jarnik, refined by Hensley [1996], we have 
f(a) = 1 — := 1 — but only if m — and therefore a — is not too 
small. 

The result is partially hinted at by Cesaratto and Piacquadio [1998], 
Piacquadio and Cesaratto [2001], and Piacquadio [2004], and in Piacquadio 
[2004] it is empirically shown to be computationally correct within relatively 
small percentage errors. 
Note: The value of 

just quoted, responsible for /(a) in the F — B case, when a is not small, 
can be refined a bit. Let us remember that the exponent "-2" comes from 



Eq. 35 the exponent is a = —2dii{Ek) — — 2 if k (and k m , and m, and a) is 
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large. Remembering also (see the end of Sec. 5.2) that - r^pjp is the generic 
F — B contractor, and that djj(Ef e ) with E^ restricted by i^Zj a j)/ n < m 
becomes f{a) in this section, we finally have 

(j + l)-2/(a) 

A, = (17) 

[ i yj/Ca) 

which can be written as C7r(j+ g , where the last £ normalizes the in- 
troduced factor (^-)-^ a ) • So we have Xj = - — contr ^ ctor ] /( } _ Notice that this 
value has much in common with the critical A's for the Euclidean case: for 
equal lengths we had 

[j th contractor] ex P° nent 
Aj - , (18) 

as for equal probabilities, and again for the F — B measure. The difference 
is in the value of the exponent: /'(a) for equal lengths in the Euclidean 
case, f'(a) for equal probabilities, same case, f(a) for the F — B one. . . We 
will return to these apparent differences later on, in Sec. 6. For now, we 



want to stress the universal character of Eq. 18 where Euclidean and F — B 
measures intersect. 

5.4 A better expression for (a, f(a)) 

We want now a more accurate expression for a and /(a) for the F — B 
measure. With the same notation as in Sec 5.3, we want to extremize 

A l •■■ A j A / 



i og {{([c2^...(k + i)^ }} ~ A ( io g (i/{{}» " «). with condition £A 3 - = 1, A, 
the variable. This is done in the Appendix to Section 5.4. The critical Xj 
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are 

A<= g±iqv^ (19) 

Now this value of Xj seems to be very different from those obtained in 
Sees. 3 and 4, and from those in Eqs. 17, 18, and those from previous 
work [Piacquadio and Cesaratto, 2001]. We will show the corresponding 
connections in Sec. 6. 

5.5 The information dimension for the F — B measure 

In this section, we will find the value of a for which f{a) = a and f'ip) = 1, 
showing that this entropic or information dimension is the universal value 
0.870 ±0.0004 found by Bak and others [Bak, 1986 and references] to be the 
approximated box dimension of the fractal underlying the Cantor staircase 
for the circle map, in frontier with Chaos. 

Equating f(a) and a we obtain -| ^§*T^ = ^ lo^i'..) ' which im " 
plies Tjji^j lQ g ^+(log 2)j\ j ) = or £\ A,(log A.+log V) = £V A, \og{X0) 
0, so, if we write Xj = i we have ^ • Xj = 1 and f(a) = a. 

1 _ (j + l)2r(A) /2 A(j-l) 



3G 



now reads 2 , 



With this particular value of Xj , Eq. 

which can be rewritten (with £ always the corresponding normalizing sum): 

1 (j + l)2r(A) 



2i 2 A i£ 



(20) 



Or ^ = 2( A -^' = const. (j + l)2(Aa-/(a)) = const (j + 1 )2/(a)(A-l) j since 

a = /(a). So we have 

2( A - 1 ^ = const(j + l) 2 ^ A - 1 \ (21) 
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with 2/(a) a number strictly between and 2. Now, if A ^ 1 we obtain 



from Eq. 21 2 J = const (j + l) 2 -f( a ), obviously an absurdum as j grows, so 
we confirm that a = f(a) implies A = 1. On the other hand, let us assume 
that A = 1 above, in Eq 



20 



We are left with = const y — '—^ , which 



implies that (j + l) 2 ( Q- /( a )) is a constant Wj, an absurdum unless a = f(a): 
so a = /(a) iff A = 1. But A = 1 and a = f(a) mean r = 0, which seems to 
be in agreement with the Euclidean cases as the condition that characterizes 
the concentration of the measure i.e. the information entropic dimension. 

For this case, in which Xj = i-, we have the corresponding f(a) = 
— j> i Al /of ^n-L'" n to be 0.87038. . ., the Hausdorff dimension of the subfrac- 

2 log(c2 A i3 A 2...) ' 

tal which concentrates the F — B measure. Notice that this number lies 
in the interval 0.870 ± 0.0004 quoted above. The more restricted interval 
0.870 ± 0.00037 [Weisstein, 2005] for the box dimension of the fractal as- 
sociated with the Circle map staircase, would differ from 0.87038 ... in one 
unit in the 5th decimal, an error that arises from the use of the simplified 
Besicovitch approximation — bound to be "very good" indeed, according to 
Good [1941] — which does not take into account the order in which the par- 
tial quotient coefficients dj appear in the cumulant q n = q n (a\ . . . , a n ), but 
only their values. 



Observations. The formula for the key Aj's shown in Eq. 36 is much 
more complex than those for the Euclidean cases. An adaptation of the 
reasoning in Sec. 4, in order to prove that A = f'(a) in the F — B case has 
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been, so far, elusive. That is why we showed in some detail that, at least 
in the case of the subfractal that concentrates the F — B measure, A does 
act as f'{a). These efforts are necessary, when we recall that the validity 
of the thermodynamical formalism has been proved only for the Euclidean 
measures [Cawley & Mauldin, 1992; Riedi & Mandelbrot, 1997; 1998] and 
only semicomputationally for the F — B measure. [Piacquadio & Cesaratto, 



There are old and new results on the (a, f{a)) spectrum of the F — B 
measure, and in the next section we harmonize the — only apparent — cor- 
responding discrepancies — not all details included, for obvious limitations 
of scope and space, some fine brushings are left to the reader. The reader 
only interested in following the thread of the argument on d ~ 0.870 may 
skip Sec. 6. 

6 Relating the Key Aj's 

We seem to have two — apparently — very different expressions for the key 
Aj's in the case of the F — B measure, which are, in turn, quite different 
from the key Aj's corresponding to the Euclidean case. Let us study these 
apparent discrepancies. 



2001] 



For the F—B case, the value of Aj from Eq. 36 is A 



)- T /2(-'- 1 ' A 



'j — 



e 
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and we want to connect this result with the value 

1 1 



(j + I) 2 Cn 

j : 1 — > k m = (e C7r ) m := e Ba , since average m of the a,j is proportional to 
a, all according to Sec 5.3 ; f(a) = 1 - := 1 - A and 5 positive 
constants, B > 1 — and let us recall that this result was valid when to, and 
therefore a, was not small. 

Let us continue to assume that A is the derivative of f(a). Then A = 

AB _ const _ j f we recall that • < k we have A • < const . gQ the yalue 

is, 



/in 



of Aj, the exponent of 2 Ajf above, is bounded, ...so Aj from Eq. 
essentially, ( (j+i)2 )~ T , normalized. 



3G 



Now, let us have a closer look at the other expression (Eq. 22 ) for the 



i 2 ( i )/(«) 
key frequency: Aj = jj^rpj2^~ = (J g 1) ? the approximant of 



(i'*' 1 contractor)^") j- -r-i 

— g according to Eq. 



35 



The exponent — r in the expression 

would be, 



36 



above, "( )~ T , normalized", is — r = f(a) — Aa, so Eq 
essentially, ( jj^[yz)^ a H (J+Tp)^ ' norma hzed; j : 1 — > k m . We want to 
analyze, therefore, the behaviour of the discrepancy between expressions 
[IT] and Q i.e. ((j^)"^ = ((j + l) Aa ) 2 , j : 1 ^ k m = e Ba and A = 



= being proportional to to. So (^w) _Aa = (j + l) 2Aa < 

(k m + l) 2Aa ^ (e Ba )^ a = e ^ const » 1. If j does not grow, still the 
exponent Aa tends to zero and, again, ( (jjip )~ Aa ~ 1. So both expressions 
of the key Ay's are very much like (^pjp) -1 ", normalized. Finally, if we 
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recall that const. , -2-,^ is the generic contractor in the F — B construction 



is, exactly, the value for the key \j in the Euclidean case. 

7 A Statistical Version of the Farey Tree 

By connecting the cases where a segment is measured with a common Eu- 
clidean ruler, or by the F — B probability \/2 N , we tried, so far, to establish 
a connection between Euclidean and F — B measures, by means of their 
corresponding multifractal analysis. The differences between the two mea- 
sures are considered to be deep and are briefly pointed at in Sec. 9. Yet, 
the thermo dynamical algorithm — the multifractal spectrum — reveals, on a 
closer look, their inner links. We propose to deepen these links. 

Let us suppose we are studying, empirically, the geometry of a fractal 
in a unit segment given by, say, a certain dynamical system, so we know 
the step N in which we are. Further, let us suppose that the fractal is 
— once constructed, as N grows — a ternary like that of Cantor, a typical 
self-similar "Euclidean" case in the sense described above. The subdivision 
of segments seems to correspond, empirically, to a left-right process, so we 
know that in step TV" we have a list of 2^ segments. Their length seems to 
diminish exponentially, like ^,^4 > 1, but we are not sure of the value of 
A. We are not so much interested in the value of A, but on that of log A, for 
we know that would be the dimension that we are trying to estimate. 
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In order to estimate log A (if we are in the ternary of Cantor, log A should 
be log 3, but we are measuring experimentally) we take all 2^ segments in 
the Nth step, we take their reciprocals (so we would have 2^ segments of 
length A N , roughly), we take their logarithms, we divide said logarithms by 
N. . . and we take the average of all these values, for as large a value of N 
as we can handle. That should give us a stable value converging to logyl, 
log 3 if we were in the ternary of Cantor. 

We propose to do such a calculation for the 2 N intervals in the Farey tree 
partition: we will take their Euclidean lengths, take their reciprocals, take 
their logarithms, divide them by N, and average all these values. This will 
be our log A, and will be the dimension of the Euclidean statistically 
self-similar version of the Farey tree. 

Let us recall that we have |^ (ai . . . a n ) = - + 1 1 — in the step a\ + 

'"an 

. . . + a n = N of the F-B partition. We estimate q n (ai ■ . ■ a n ) according to 



Eq. 12 as [c2 Al 3 A2 . . .] n , where Xj is the proportion or frequency in which 
a coefficient equals j. Therefore logq n = n[logc + Ailog2 + A2log3 + 
. . .] = nlogc + Zilog2 + ?2log3 + where lj is the total number of co- 
efficients a,k = j. Then logq n = nlogc + ^jl°s( a j + !)• We are in 
step N = a\ + . . . + a n . We also recall that we estimated length of seg- 
ments as \, so, if we take the reciprocals and take logarithms we obtain 
2 log q n = 2{n log c + ^ ■ log(aj + 1)}; Y^j a-j = N. Before dividing by N, we 

log A 



will take averages of these values, in order to obtain log A, ^4 the dimen- 
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sion of the Euclidean version of the Farey partition: we have to average the 
index "n" in a certain iV-step; in order to average the values ^ ■ log(aj + 1) 
we have to count first how many coefficients a,j = 1 we have in step N, how 
many cij = 2. . . until aj = N, which happens only once in that step. Then 
we can take averages and calculate log A. The whole counting-and-averaging 
process, long and tedious, is done in the Appendix to Section 7. The result 
is lo gj 4 = log C +¥ + ... + i^il) + ... . 

8 log 2 / log A is the Information Dimension for the 
F — B Measure 

We want to compare this "Euclidean" dimensional version log 2/ log A of the 
F — B measure with the entropy or information dimension for the F — B 
measure in Sec. 5.5: f{a) = (- \) logc+Ai i^jf^fe^ log4+ .„ , for \ 3 = ±. 
The denominators coincide. We have to compare (— |)(Ai log Ai + •••); Xj = 

1, with log 2. We have log A-,- = log i = — jlog2. So — |(AilogAi + 
...) = i(!l og 2 + i J 21og2 + ^31og2 + ...) = ^{i + ^ + J^ + ...}. 
Let us calculate the expression within brackets. The Taylor expansion of 
jj^xyi is 1 + 2x + 3x 2 + ... + (&+ l)x k + — For x = \ we have then 
1 + \ + §i + ■ ■ ■ + + • • • = ( 1 A)2 = 4 which implies, dividing by 

2, ^ + ^ + + + ^ + = 2; so 2 is the value of the expression 
within brackets, and log 2 is the value of the numerator of f(a) above, which 
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means that the expressions for f{a) and above coincide, so 0.87038 . . . 
is both the information or entropic dimension of the F — B measure and the 
dimension of the Euclidean version of the Farey tree partition. 

9 Geometrical considerations, Conclusions and Con- 
jectures 

9.1 Geometrical considerations 
9.1.1 Introduction 

For the content of this section we refer the reader to The Geometry of Farey 
Staircases [Piacquadio, 2004] and to the corresponding references quoted 
there. 

There is a one-to-one connection between F — B in II and a certain non 
Euclidean geometry. Though we work on /, the F — B interpolation is valid 
in any interval [n, n + 1] , n G Z. 

Let H = {z = x + iy/(x,y) G H 2 ,y > 0} be the upper half plane. We 
draw in H the upper half circles (centre in R) with endpoints in a pair of 
adjacent rationals |, y in any Nth F — B partition. That is, we trace upper 
half circle (centre ^) joining and 1, then arc joining with |, then | 
with 1, . . . etc. in the Nth partition we trace 2^ small arcs joining adjacent 
rationals as endpoints. These arcs are geodesies in H. The three geodesies 
joining ^ with ^ (adjacent in an Nth F - B), ^ with |±^, and with 
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y (in (N+l)th F — B), form a triangle in H. We have an infinite number of 
such triangles and, in the so-called Hyperbolic area measure, they all have 
the same area. 

A rigid hyperbolic movement in H is, by definition, a transformation 
pf±§,£ S H, a, a', b, b' in Z, det 



z 



a' a 



= 1. The set U of these 



movements can be seen as the multiplicative group of 2x2 matrices with 
integer entries and unit determinant. The triangles above, do not only 
have the same hyperbolic area, but are transformed into each other by rigid 
hyperbolic movements: by elements in U: they are — hyperbolically — the 
same triangle, moved here and there, to and fro. We do likewise in any 
interval [n,n+l],n£Z. 

To the arcs described above, let us add vertical lines (n, oo) with endpoint 
nEZ — which are also geodesies in H, the centre of the circle at infinity of 
R. On top of unit arc joining and 1 — we will call it unit arc hereinafter — 
we have now another triangle, the sides being vertical line (0, oo), unit arc, 
and vertical line (l,oo), vertices being oo, and 1. The same happens on 
top of arcs joining n & n + 1, n in Z. These new triangles have the same 
area as those above, and are interchangeable with them by elements in U. 
All these non-overlapping triangles — with finite or infinite vertices — cover 
H: they are a regular tiling of H, and we will call it T (T for triangle and 
T for tiling). 

In Sec. 2 we saw that, if i = \a±, . . . , a n , . . .1 and \ai . . . a n ] := then 
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length of segment T 211 , 2™+Il j s — I — which implies det 



and 



Pn+1 Pn 



£ f/. This means that 2:1 and 211:111 are adjacent fractions 



in some F—B partition, for all adjacent rationals | , y- in all F— S partitions, 



| < ^,have 







( a > 


\ 


a' a 






a 




= 1, i.e. 






b' b 









G U: there is a common structure 



in charge of F — B, continued fractions, and rigid movements in Hyperbolic 
Geometry; the algebraic group U being the common underlying principle. 

( , \ ( ■ ^ 



if 



a a 



b' b 



€ U, < ¥ < ?r < 1, then 



a a 



b' b 



applied to unit 



segment shrinks I into [|, ], yielding the F — B interpolation between 
adjacent | and p-; the length of the shrunk interval is - Second row 

entries b and b' are non-zero and positive. Ditto when working in [n, n + 
1] instead of /, n G IN. When n is negative, such entries are non-zero 
and negative. But other elements in U can have b and b' of different signs 



or zero, e.g. u 



n + l -1 
1 



and u = 


( , 

1 — n 


\ 

n 


) 


v - 1 


1 j 



,n € Z. In 



both cases the denominator 6(6' + 6) above is zero: 0(1 + 0) and 1(— 1 + 
1), respectively. Element u* transforms unit arc into vertical line (oo,n), 
whereas u transforms unit arc into vertical (n, oo) — same line with different 
orientation, so u and u* mirror each other — and unit segment into horizontal 
[n, oo): so b(b' + b) = for u and u* shows that length b ^ b } +b ^ of those lines 
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— a(unit arch) and u*(unit arch) — is infinity. Every u in [/has a u* mirror, 

related to u in a technical way beyond the scope of this paper. An analogous 

1 n 



analysis can be done to the translations z — > z + n, i.e. 



[/, n € Z. So, the correspondence between F — B and [/ goes beyond 
the corresponding to non-zero-equal-sign-denominators of fractions in /, but 
extends to semicircular arcs with F — B adjacent endpoints throughout R, 
and to vertical lines with endpoints in Z, i.e. to all geodesies delimiting all 
triangles in T. Notice that every such geodesic is obtained by applying each 
element of U to unit arc — which is the rationale for restricting the work to 
unit interval in the next sections. The u/u* mirror ambiguity is avoided by 
joining smaller with larger values: to 1 in unit arc or segment, n to oo in 
the infinite lines. Other regular F — B tesselations of H aim to take care of 
this apparent ambiguity, but we stick to T, simpler to work with, and which 
embodies all geometric and metric properties of H, as well as defining, via 
the endpoints, the F — B partition on H; which is the reason we have used 
the terms F — B measure and hyperbolic measure in R as interchangeable 
in the literature. 

9.1.2 Equal F — B measure of all intervals in Nth F — B partition 

Let i = [01,02, • • ■] £ I- Two matrices L and R — for left and right — in 
U can be constructed such that, applying L a\ times to I means to F — B 
interpolate I a\ times, each time choosing the left interval in order to F — B 
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interpolate further. Ditto for R and 02... and so on. So [01,02,...] = 
lim n ^ 00 L ai R a2 . . . L 2n+l I. Now L and R — left and right — are like, say, 
vectors (0, 1) and (1,0) in R 2 — horizontal and vertical: they carry the same 
weight, have the same "right to be present". So, e.g. L ai R a ' 2 L a3 will have 
the same F — B or hyperbolic measure as L a i+ a 2+ a * or LR ai+a2+a3 ~ 1 : all 
words written with N letters L and R have the name weight or probability 
measure f° r eacn of the 2^ intervals in step TV of the F — B partition. 

9.1.3 A deeper connection between H and F — B 

Let us consider a vertical geodesic Gi in H with i = [a\ . . .] € I as end- 
point. It cuts an infinity of triangles in T. Let us trace Gi with a finger 
at its left side, from top to bottom. When crossing a triangle through a 
thin part (only one vertex at left of Gj) we write T for thin, otherwise we 
write F (for fat) — the tile at the very top of Gi is T, for technical reasons 
beyond this paper. We obtain an infinite word, letters T and F: T ai F a2 . . . 
identical with L ai R a2 ... in last section. This fact tightens the connection 
between continued fractions, F — B, and H. The main point here is that 
by naked eye observation, tracing Gj with a finger, we can write directly 
any i in its continued fraction, without any calculation. Let us recall (Sec. 
5) that irrationals with the same a-concentration are those with, roughly, 
the same average over the Oj values: this can be verified by looking at Gj's: 
cardinality of tiles in T with adjacent T's or F's should be — statistically — 
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the same. Also, knowledge on the a/s of an irrational i, implies knowing the 
classification of said i (Bruno, Jarnik, Liouville. . . ) needed by physicists to 
study circle maps or optoelectronic phenomena [Piacquadio & Rosen, 2007]. 

Now: suppose that we have the ordinary half plane R 2 with an ordinary 
regular tiling, all tiles interchangeable by rigid Euclidean motions. Notice 
that no geodesies — vertical or otherwise — crossing the tiles with endpoint 
in an irrational t £ E will yield those tools to classify said % according to 
the criteria needed by physicists, whereas any geodesic in H with endpoint 
in i — not only the vertical Gi — will yield such classification. 

9.1.4 Fundamental differences between Euclidean R 2 and H 

So far, we stressed the tight connection between F — B, continued frac- 
tions and H (cum U cum T), with an emphasis in H. And, at the end 
of last section, we pointed out like a divorce between upper half planes 
R 2 and H. Such differences run deep indeed: we can have R 2 regularly 
tesselated by triangles, squares, hexagons. . . period, whereas it is a most 
satisfying exercise to transform H into the Poincare circle P, to choose, say, 
five or eight consecutive geodesies, and tesselate P (hence H) with regu- 
lar pentagons, octagons,. .. etc — an impossible endeavor in R 2 . Opposite 
characteristics are easy to observe even at the level of R: when we write 
x = 0, aq . . . a n 00 . . . € /, we know that x is rational, whereas in F — B 
hyperbolic style a rational in / is written [ai, . . . a n , oo, oo, . . .]. Likewise 
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x = 0, a±a2abab ... is rational, whereas [a±, a,2, a, b, a, b, . . .] belongs to the 
most irrational category. . . 

9.1.5 Analogies between the two measures, Euclidean and Hy- 
perbolic 

The list above of apparently irreconciliable differences between the two mea- 
sures is by no means complete, for many more are pointed out in the lit- 
erature. Some analogies, instead, have been noticed in [Piacquadio and 
Cesaratto, 2001], and they begin to appear, obscurely, through multifractal 
analysis. 

In Sec. 7 we take the Euclidean length of the 2 N F — B intervals in step 
N with a common Euclidean ruler. We obtain a list of 2^ values, of which we 
take logarithms. Some values are larger, some are smaller, so we take their 
average, which yields a single statistical contractive value \ £ (0, \), so 
is the statistical self similar dimension of the F — B partition. In the ternary 
of Cantor \ = \, and we have 2 N subfractals — in segments of Euclidean 
length — interchangeable by rigid Euclidean movements, ... all of which 
happens with the single contractor ^ above: it yields a — Euclideanly — 
self-similar fractal, a statistical counterpart of the F — B partition. But in 
Sec. 8 we learn that its dimension |°g ^ is the same d ~ 0.870 . . .: here is a 
deep contact between Euclidean and Hyperbolic geometries. 
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9.2 Conclusions and Conjectures 

The two measures, Euclidean and Hyperbolic meet in a very specific di- 
mension: 0.8700 ± 0.0004. This value of f(a), the entropic or information 
dimension, corresponds to the fractal where the Hyperbolic measure is con- 
centrated — whereas the dimensional Euclidean perspective "sees" the Farey 
Tree partition as having this specific dimension, instead of dimension 1. This 
universal number, therefore, is strongly perched on, and comfortably accom- 
modated in, the intersection of the two measures. How does it appear in the 
dynamics of the Circle Map? For just a moment let us suppose we under- 
stand that the Circle Map acts as a black box: the input is the "y" -vertical 
axis in the Devil's Staircase associated with the map: the entire unit seg- 
ment is there, the input is, dimensionally, 1. The output is the selected 
subfractal O, C / in the horizontal axis (associated with the circle map stair- 
case) of dimension 0.870 .... This black box seems to act as a dimensional 
spaghetti percolator: the output, what is retained, is, Hyperbolically, that 
set where the measure is dimensionally concentrated, yielding full informa- 
tion on such measure. This would be the thick fat spaghetti, whereas what is 
lost, what oozed through the percolator holes is the very small stuff: herbs, 
salt, fine flour, seasoning, small particles that came with the spaghetti in 
the input,. . .which do not yield much information, do not concentrate the 
measure. 

From the Euclidean point of view, the whole of the input is dimension- 
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ally retained in the percolator, for the Euclideanly self-similar version of the 
Farey partition has exactly this 0.870 . . . dimension. 

Let us assume we accept that the circle-map-Devil-Staircase black box 
acts as such a percolator: it retains the concentration of information. Then, 
the universal character of this numerical constant might be clear: changing 
the "sine" function in the circle map by another reasonably smooth func- 
tion that draws the circle, would mean changing a percolator by another of 
a slightly different form, say, an enamelled one with little circular holes, by 
a wire net one with adjacent square holes: the same spaghetti would remain 
trapped, the same output would be obtained, the same tiny particles lost. 

Why and how the circle-map-Devil-staircase black box acts as such a 
measure percolator, however, still remains, for us, a mystery. 
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11 Appendix to Section 3.2 

We have to extremize -^^(Ai log Ai + . . .) - Aj^^(Ai logpi + . . . - a) + 
fj,(X\ + ... — 1), so the derivative with Xj as variable: — i Q g no (1 + log Xj) — 
A E^b log ^' + ^ = i m P ne s 1 + log Xj - A log p.,- + fi* = 0. 

Subtracting the equation corresponding to Ai we have log ^ — A log ^ = 

°> or it = (^") A ' J = 2 ~* n o- Tnen Aj = ^kpf, j = 1 -> no, and from 
2J • Aj = 1 we obtain Ai = * A = hence 

X - Xl n A - P l- 1 t, a - ^ Vi f23^ 

Aj - ^ - T ^ " £ Vj (23) 

The condition — log 1 Wo Ylj Aj logpj = a permits knowing the value of A: 
A = A(a), or a = a(A) (24) 



With 



/(a) = i Y] Aj log Xj) a = V Xj logpj; Xj = -jr (25) 

log n ^ log n ^r 1 £ 



3 3 



we can calculate f'(a) = 4^ = 4^44 = -Frt r J " > J = g J l/ J since 
Aj- = 0. Therefore f'(a) = ~^^^^t^- = A, again since £V A^ = 0. 
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12 Appendix 1 to Section 4 

We have to extremize the function Al ~ A ( \ a)+/u(Ai + . . .- 

"1 ^5 Ax iOg Cl~r"... 

1). The derivative of this function, A-,- as variable, equated to zero yields: 

(l+logA j )(Ailoge 1 +...)-(A 1 logA 1 + ...)logc J . (-logn ) , , .. _ n 

(A 1 logc 1 + ...) 2 + A (A 1 logci+...) 2 ° gC ? + ~ U ' nence 

log Aj(Ai log ci + ...)- (Ai log Ai + . . .) log Cj - A log n log cj + fj,(Xi log ci + 
. . .) 2 + (Ailogci + ...) = or logAj(Ailogci + . . .) - log Cj(Xx log A x + . . .+ 
A log no) = a constant independent of j. Subtracting the corresponding 
equality for j = 1, writing Ai log c\+. . . = (), and Ai log Ai+. . .+Alog no = [] 
for short, we have () log & - log g- [] = 0, or & = (g)0/O ; Aj = -fecJ /0 
with j = 1 — ► no- 

As in Sec. 3 we use Yj Xj = 1, we obtain Ai and then Xj, with the result 

J/0 r s 



which bears a resemblance to Eq. [25| 

Next, we want to calculate /'(a) = f-. We write /(a) - AllogAl 



da' J \ J Ailogci+... 

S for short. Then a = =%*L. Therefore /('a) = £ = g§ = 

{num'dcn— num- dcn'j/dcn 2 _ 1 [(l+logAi)A / 1 +...](Ailogci...)-(AilogAi+...)(A / 1 logci+...) 
- log n ( - 1 /den 2 ) den' logno Ajlogci+... 

Since £ A;- = Owe have /'(a) = ^{^^(Ax log c 1+ . . .)-(Ax log A 1+ 
...)}. Now, log Aj = S log Cj - log £ implies A^ log Ai + . . . = Y^j s log CjX'j - 
log £ Y,j Xj = S (X log ci + . . .) . Hence, 

/'(a) = r^{H(Ai log Cl + ...)- (Ai log Ai + ...)} = 
log n 

{^(Ailogd + ...)- (AiHlogd + . . .) + log£(Ai + ...)} = ^-(27) 

log n log n 
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13 Appendix 2 to Section 4 

Variable a, the natural independent variable of / is 

1 , s 

a = -, (28 

Q 

a the "old" concentration from Sec 3. The inversion formula of Riedi and 
Mandelbrot [1997 and 1998] says that the new inverted spectrum / is related 
to the old / : /(a) = a/(±). From this we have /'(a) = /(I) + a /'(i)^ = 

/(a) ~ a/'(a)' and from above this becomes 3 = /'(«) = /'(£) = /H ~ 
af'(a) = f(a) — aq = —(qa — f(a)) = —r(q), that is 

- r{q) = q Vg (29) 

where r and g are "old" parameters. Applying again the same criterion we 
have (g) = —f(q) = — f(— r(g)), that is q = — f(— r(g)). But f = f implies 
q = q, so we have 

-f(-r(q)) = q Vg or - r(-f (g)) = q Vg (30) 

Let us go back to Eq. [27} in our new notation: 

df log£ -Nlog(cf + ...) -log(cf + ...) N 



da log n ^log^ log ^ 



lo g Eri+.-=JV C ? ri • • • c n r n ° log Er 1 + ...=iv( c ? ■ ■ ■ c «o° f ,„x , Q1 v 

= ( 31 ) 



log(i)* log(length) 
for we have the case "equal lengths and different probabilities" given by 
contractors Cj, . . . so we are in Sec. 3 with the very definition of t. 

Now, let us focus on variable H = Jl = io g y ...+Aiogn = A ^ ogAl +- - 

' () Ailogci+... AilogCi+... 



14 APPENDIX TO SECTION 5.2 41 



A Al |l° g r r° = /(«) - Aa = -(Aa - /(a)) = -r (A), so Eq. 



31 



now reads 



— t(3) = — t(— f (A)) = A VA from Eq. 30 Hence, we have our derivative 
f'(ct) = A, as was the case of "equal lengths" in Sec. 3. It means that 
q = A in Sec. 4, as q = A in Sec. 3, in both cases A being the Lagrange 
indeterminate coefficient linking f(a) with a. 

14 Appendix to Section 5.2 

We need, now, to estimate constant c G (0,1). We will adapt a reasoning 
that we used elsewhere [Piacquadio, 2004] in order to apply the methods in 
Sec. 4. 

Let Ek = {i = [a\ . . . aj . . -]/aj < k V?}, k £ IN. We will estimate 
the Hausdorff dimension dniEk) by considering finite sequences [ai . . . a n ], 
aj < k (later n will tend to infinity), and considering — as above — the Xj 
as the frequency in which the a's are equal to j. For a certain choice of 

log(l/A Al A Afe ) n 

Ai . . . Afc we have the corresponding dimension given by 1 — k — 



log(c2 A i...(fc+l) A fc) 2 ™ ' 

simply repeating the processes above. So, dn{Ek) will be obtained by finding 
the key set of frequencies Aj for which 

_1 Ai log Ai + . . . 1 num 

21ogc + Ailog2 + ... + \ k \og{k + 1) "~ 2 den ^ ' 

reaches its maximum. 

Again ^(^j^) = implies (1 + logAj)den = numlog(j + 1), that is 
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1 + log Xj = ^ log(j + 1), or, from Eq. [32 



1 + log Xj = -2d H {E k ) log(j + 1) := alog(j + 1). (33) 

Following already well trodden steps, we obtain log ^ = log[ ^ 1 ^ ] a from 
above, from which Xj = ^t(j + l) a . With ^ - Xj = 1 we obtain 

2 a (i 4-l) a 

Ai = -, and Xj = u - £ ' . (34) 

Now let us examine the value of 

a = -2d H {E k ). (35) 

As k —* do, d}{{Ek) must tend to unity, as E k tends to encompass every 
i = [ai . . . a,j . . .] regardless of the size of the a's. In fact, Jarnik [1928; 
1929] proved 1 — dniEk) = 0(f)- Therefore, as k grows, a = —2. Hence, 
£ ZTU + !)~ 2 = Thus Eq. 



34 



becomes Xj = 2 , 



k large, and Eq. 32 becomes (always k large) —2dH(E^) = —2 = 



iogc+^E,(i+i)- 2 iog(i+i) - " iogc+EEE trom whlch " 2 log c 



£ SO' + 1) 2 lo g c tt = - log c n , that is c 
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C-TT . 



With the same notation as in Sec 5.3, we want to extremize ~* 

lo, 



Iog( , 1 A ) n 

\ A 1 \ k ' 

log{{([c2 A l...(fc+l) A *=]»)2}} 



^( log(l/{{}>) ~~ a )' w ^b condition ^ Aj = 1. Proceeding as in the Euclidean 
case, we have to find extremes of — \ ^ 1 ? g r,v; + T — A '° S 2 .j?/^ ■ that is, the 

' 2 log(c2 A i...) -21og(c2 A i...) ' ' 

extremes of Al tog^-f^ := • ^ e e< 4 ua te the derivative of this 
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function (variable Xj) to zero, and the difficulties in the (apparent) differ- 
ences with the Euclidean case begin to appear: ((log Xj + 1) + j'A log 2) den = 



numlog(j + l), so 1+log A j +jAlog2 = ^log(j + l), with ^ = -2/(a) + 
A raSfc = -2/(«) + 2Aa = -2(/(a) - Aa) = 2r(A). Therefore, 



proceeding as before, we have log ^ + (j — l)Alog2 = 2r(A) log ^ ' , or 



log ^ = log(i±i)MA)_l og 2 A 0- 1 ) which means Xj = Ai(i±i) 2r - , and 



the equality is valid for j = f as well. With J2j = 1 we obtain, as before, 
the value of Ai and then that of A,: 



16 Appendix to Section 7 

In order to calculate log A we need to closely study the nature of the re's 
and Oj's in a certain N step. 

16.1 Rewriting the tree 

We start with the first interpolations of the tree: 



(j + l)2r(A)/ 2 A(i-l) 



(36) 



£ 



o 4 
l ' l 



i 

2 




1 

3 



2 
3 



1 

4 



2 

5 



3 
5 



3 
1 



Let us express, in terms of continued fractions, the values |, |, | and |: 
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| = t = TZT' so the a fe involved are 1 and 2; % = \ = -j~r, so the 

2 1+ 2 2 2 +2 

involved are 2 and 2; | = \ = = = - — , so afc are i j i j anc [ 

3 1+3 l+ f 



3 

rewritten thus 



2. And | = 1 = TXT, with a& being 1 and 3. So the tree above can be 

[2] 

[3] [1,2] 

[4] [2,2] [1,1,2] [1,3] 

So the step can be seen as the sum of the a& involved inside brackets in each 

horizontal line: 

N=2 [2] 



N=3 ... [3] ; [1,2] 1+2=3 

N=4 [4] ; [2,2] ; [1,1,2] ; [1,3] ...2+2=1+1+2=1+3=4 

and so on. Let us examine the minitree 




We observe that the cypher 1 in [1,2] appears in both the daughter branches, 
whereas the last a& in [1, 2], i.e. 2, unfolded thus: [...2 — 1,2] or [...1,2]. 

This is general, as we see by examining the other minitrees, or by extend- 
ing the tree to N = 5, 6, The process: [a±, . . . , a n _2, a n -i] m step N—l = 

oi + ... + a n _i generates [ai, . . . a n _ 2 , a„_i + 1] and [a±, . . . a„_ 2 , a n _i - 1, 2] 
in step N. 
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We need to estimate the average of all lo |^" in a certain step N. The 
tree above has only the new elements which appear in step N; we will call 
this the restricted tree — restricted only to the new elements in the step. 
The tree with all elements will be called the complete tree. The elements in 
the N th horizontal line of the restricted tree will be the restricted elements 
in step N. 

16.2 Averaging index "n" in an iV-step 

Since log q n is estimated by n log c + J2j l°g( a j + 1) we will start by averag- 
ing all nl ^ gc involved in step N. We start by adding up all the "re's" in a 

restricted iV-step. Let us enlarge the restricted tree a bit more: 
N = 2 — >[2] 

AT = 3— [3]; [1,2] 

AT = 4— [4]; [2,2]; [1,3]; [1,1,2] 

N = 5— [5]; [3,2]; [2,3]; [2,1,2]; [1,4]; [1,2,2]; [1, 1,3]; [1, 1, 1,2] 

The corresponding values n(N), for step N, i.e. the lengths of integers 
afc inside brackets are: 
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N n{N) 

2 ... 1 

3 ... 1,2 

4 ... 1,2,2,3 

5 ... 1,2,2,3,2,3,3,4 

Rearranged, these numbers are: 

N n{N) 

2 ... 1 

3 ... 1,2 

4 ... 1,2,2,3 

5 ... 1,2,2,2,3,3,3,4 

6 ... 1,2,2,2,2,3,3,3,3,3,3,4,4,4,4,5 

By simple observation we see that the lengths n(N) vary from 1 to N — 1, 
and that their repetition follows the combinatorial numbers in the Pascal 
triangle of order N — 2. So the sum of all the n(N) in the restricted N step 

is ( V) • 1 + (V) • 2 + . • • + (S3 (iV - 1) = Ef=o 2 (V) (J + 1), and we 
claim that this sum is N2 N ~ 3 . It is obviously true for N = 2, the first N of 
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our account. Let us assume the validity of 

J:( N : 2 )U + 1) = N2^ (37) 

3=0 V 3 7 

for a certain value of N, and let us infer the validity of said equation for 



N + 1. For short, we write N — 2 = k, so Eq. 37 becomes Y^j=o (j) + 1) 



(fc + 2)2 fc ~ 1 . We have 

fc+l /, , 1N k 



E (* + 1 )0- + i) = i + ^('= + 1 )0 + i) + (fc + 2) 

7=0 \ J ' j=\ \ J / 

i + dGK-> +i)+( * +?) 



A- 



•+eO' +1)+ £G-i> + eG-i) +( * +2)= 

i=i VJ/ j=i v 7 i=i v 7 

E C) w + » + '~E 1 G * > -i) + i] + (*+i)+E (*) - 

j=o VJ/ .7-1=0 v 7 i=o VJ/ 

2 ^ Q (j + 1) + 2 fe = 2(jfc + 2)2 fc ~ 1 + 2 fc = (Jfe + 3)2* 



which is Eq. 37 with k replaced by k + 1. Therefore, the sum of all n(N) in 
the restricted step N is N2 N ~ 3 . 

Next, we need the sum of all n(N) in step N: 

2 • 2 2 - 3 + 3 • 2 3 ~ 3 + . . . + N2 N - 3 = ^{2x 2 - 1 + 3 • x 3 ' 1 + ... + Nx N ~ 1 }\ x=2 
= \{x 2 + x 3 + ...+ x N )'\ x=2 = \{x 2 {l + ... + x N - 2 )}'\ x=2 = 

Now, we are interested in computing the average, in step N, of values n ^ : 
we need to divide j2 N (N— 1) by iV and by the total number of tree elements 
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[5] ffl [2?] [2,1,2] [M] 

/\ /\ /\ /\ /\ 



n 



H.H^P^PWlPflP^lP^pwaMP.A?! PW] [1^,1^1114.4] [i,uj][i,i,io][i.i,i,U] 



etc 



Figure 1: Diagram D\ 

in step N, which is 1 + 2 + . . . + 2^-2 = 1 = 2 JV_1 - 1. That is, the 

average we look for is \2 N (N- 1)/N(2 N ~ 1 - 1) = ±(1 - ^) 2 ^_ 1 ; the limit 
value when N grows is \. 



16.3 Averaging YTj=i l°g( a j + 1) i n step N 

We need now to repeat the process with magnitude Sj=il°g( a j + 1) for 
all tree elements in step N; i.e. we need to sum all log(aj + 1) for each 
coefficient aj which appears in each tree element in step N, and then divide 
the sum by N(2 N ~ 1 - 1). 

We start by considering, again, the restricted tree: see Fig. 1 (Diagram D\). 
We notice that the sum of the coefficients aj inside a pair of brackets equals, 
exactly, the value of N in which this tree element is located: the last one of 
the last row: [1, 1, 1, 1,2] fulfills 1 + 1 + 1 + 1 + 2 = 6 = N. We also recall 
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the law shown in Fig. 2 (Diagram D-i). 

[ai, a 2 , . . . , a n _i] in row N — 1, with a j = N ~ 1> 



[oi, 02, . . . , a„_i + 1] [ai,a2,...,a n _i — 1, 2] in row iV. 

Figure 2: Diagram £>2 

We start by counting the number of coefficients = 1 in step N. We observe 
that, when in a certain step a coefficient "1" appears, then it appears twice 
in the following step: e.g. ^ ^ in steps N = 4 and iV = 5 in 

[1,4] [1,2,2] 

the diagram D\. Next, we observe that, when a coefficient a& = 2 appears 
as the last one in a tree element in a certain step, then it yields an at = 1 
in the following step: e.g. ^ 2] i n steps N = 5 and iV = 6 

[3,3] [3,1,2] 

in the same diagram. Therefore, the number of a& = 1 in step N is: the 
double of the number of = 1 in step N — 1, plus the number of last 2's 
in step N — 1. But this last number, observing diagram D2 inside Di, is 
the number of elements in row N — 2, i.e. 2 Ar_4 . In a notation that, we 
trust, is natural, we can write: ajv(= 1) = 2a^^i{= 1) + 2 N ~ A . Starting 
from N = 3, where cijv(= 1) = 1, we would have: a±(= 1) = 2 • 1 + 1. In 
N = 2 there is one = 2 which is the last coefficient (as well as the only 
one). So d 4 (= 1) = 2a 3 (= 1) + 2 4 ~ 4 = 2 + 2 4 " 4 . Reiterating this law we 
have o 5 (= 1) = 2 2 + 2 5 ~ 4 plus 2 5 ~ 4 = 2 2 + 2 • 2 5 " 4 , a 6 (= 1) = 2 3 + 3 • 2 6 ~ 4 . . . 
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and ajy(= 1) = 2 N ~ 3 + (N — 3)2 N ~ 4 which satisfies the counting of l's in 
diagram D\ starting from N = 3. 

In order to obtain the total number of l's in non-restricted iV-step we 
have to add 2 Ar - 3 +(iV-3)2 Ar - 4 from TV = 3 to N = N. We have 2 N ~ 3 +(N- 
3)2 N - A = 2-2 Ar - 4 + (iV-3)2^ 4 = (2 + N-3)2 N ~ 4 = (N-1)2 N ~ 4 = \(N- 
l)2 N ~ 2 = \(N - l)^" 2 ^ = \{x N - 1 )'\ x=2 , so the sum is ^s 3 " 1 + . . . + 

X ) \ x =2 - 4(X [L + ... + X \) \ x =2 — l\X x _ x ) \ x=2 — {1\ - 2)2 

So o>n(= 1) = (N — 2)2 N ~ 3 ; an equality valid from N = 2 onwards. Here 
<27v(= 1) represents the number of l's in the non-restricted A^-step. 

From diagram D2 inside D\ we understand that it is a different problem 
to count the ajv(= 2). Again we start with restricted step N. 

A) From diagram D2 : r 1 
' 6 [ax... a n _ij 

[ai . . . a„_i + 1] [ai...a n _i - 1,2] 

we observe: 

1. each tree term [a\ . . . a n _i] introduces a number 2 at the end of 
[ai . . .a n _i -1,2]; and 

2. from diagram Di we see that no number a n -i+l in [ai . . . a n _i+l] 
is 2. So in step we have a number of "ending 2's" equal to the 
total of tree elements in step N — 1, i.e. 2 JV ~ 3 . 

B) Let us consider the case [61 2, 2]. Following the evolution that 

produces tree element [1, 1, 2, 2] in D± we observe :[1, 1, 2] — > [1,1,3] — > 
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[1, 1, 2, 2], so, case [&i . . . bk, 2, 2] in step N comes from the "ending 2" 
two steps above. . . which in turn come from the total of tree elements 
another step above: 2 N ~ 5 . 

C) The other 2's come from duplicating those in the step above: 

[ai, . . . , 2, . . . , a n _i] 

[ax ... ,2, ... , a„_i + 1] [ai, . . . , 2, . . . , a n _i - 1, 2] 

With a procedure similar to the one for counting l's this number is 

2 N ~ 5 (N-3). 

So the number ajv(= 2) of 2's in the N step of the restricted tree is 
2 jv-3 + 2 N-5 _|_ (jv — 3)2 Ar - 5 = 2 Ar ~ 5 (iV - 2) + 2 2 • 2 JV " 5 = 2 N ~ 5 (N + 2). 
This formula works from N = 4 onwards. The total number of 2's in N = 2 
and N = 3 is 2. So we need to sum 2 Ar ~ 5 (iV + 2) from N = 4 to iV = N 
and add 2 to this sum. We add 2 N ~ 5 (N + 2) in exactly the same way in 
which we added 2 N ~ A (N — 1), and we finally obtain that ct/v(= 2) in the 
non-restricted A^-step is (N + 1)2 N ~ A . 

To count l's and 2's was a special problem, but ending 3's in step N 
are produced by ending 2's in step N — 1, in a natural way, and observ- 
ing their evolution — following rules already laid out — in diagram D2 in- 
side Di, we have ajv(= 3) = N2 N ~ 5 in non-restricted step N; similarly 
a N (= 4) = (N - 1)2 N ~ 6 ; a N (= 5) = (TV - 2)2 N - J ... and so on. 

How do these quantities agree with a concrete finite row in step iV? 

Integer 1 appears (N — 2) N ~ 3 times; 2 appears (N + 1)2 N ~ A times, 
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whereas 3 appears N2 N ~ 5 times, 4 does (N — 1)2 N ~ 6 times . . .and k ap- 
pears (TV + 3 - k)2 N ~ k - 2 times if k > 3. But we know that k = N — 1 
appears only twice in row N, and k = N only once. How does (N + 
3 — k)2 N ~ k ~ 2 agree with these two quantities? For k = N — 1 we obtain 
(TV + 3 - (TV - l))2 Ar -( JV - 1 )- 2 = 4 • 2 1 - 2 = § = 2 exactly, whereas, for 
k = N the formula yields (N + 3 - N)2 N - N ~ 2 = | ^ 1: we are short by 
0.25, which does not affect our final result — it will be negligible when we 
divide the total sum by N(2 N - 1 -1). With some care we have to sum now: 
(N - 2)2 N ' 3 log(l + 1) + (TV + 1)2 N ~* log(2 + 1) + N2 N ~ 5 log(3 + 1) + (TV - 
1)2 N ~ 6 log(4 + 1) + . . . + (TV + 3 - k)2 N -^ +k ^ \og(k + 1) + . . . 5 • 2° log([TV - 
2] + 1) + 4 • 2- 1 log ([TV - 1] + 1) + 3 • 2~ 2 log(TV + 1) and divide this sum by 
TV and by 2 N ~ 1 . 

Let us add all except first and second terms above, rewriting the ele- 
ments: [TV - (3 - 3)]2 w -( 3 + 2 ) log(3 + 1) + [TV - (4 - 3)]2 JV -( 4 + 2 ) log(4 + 1) + 
. . . + [TV-(fe-3)]2 JV -( fc + 2 ) log(fc + l) + . . . + [N-(N-3)]2 N ~( N +V log(TV + l). 
Now, we divide by TV and by 2 N ~ 1 , and we obtain: 2~( 3+1 ) log (3 + 1) + (1 - 
^)2- 5 log5 + . . . + (l-^)2-( fe+1 ) log(fc+l) + . . . + (1-^3)2'( JV+1 ) log(TV+ 
1) = 2- 4 log4 + 2- 5 log5 + ... + 2-( fc+1 )log(A: + l) + ... + 2-( Af+1 )log(TV + l)- 
^{2- 5 log5 + . . . + (fc-3)2-( fc+1 ) log(fc+l) + . . . + (TV-3)2-( w+1 ) log(TV+l)}. 
Let us consider the expression between brackets: ^2 n A n , where A n = 
(n - 3)2-(^)log(n + 1). We have ^ = §^gS§ - \ if 
n — > oo; so the expression in brackets is bounded as TV grows, so when 
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we multiply it by and let N — > oo it vanishes. We are left with 2~ 4 log4 + 
. . . + 2~( 7V+1 ) log(iV + 1) as N grows. Now, we had left aside two terms: 
(N - 2)2 Ar ~ 3 log2 and (N + l)2 iV ~ 4 log3, which we have to sum, and then 
divide by TV and by 2 N ~ l . We obtain (1 - ^)2^ 2 log2 + (1 + ^)2- 3 log3 
which, as N grows, tends to 2~ 2 log2 + 2~ 3 log3. Finally, we are left 
with 2~ 2 log2 + 2~ 3 log3 + . . . + 2~( JV+1 ) log (AT" + 1) as N grows, which 

Adding all up we have the value in the denominator of : 2{^logc + 
\{ X ^ + ... + + •••}} as log A, which means logc + ^ + . . . + 

los ^ fc fc +1 ^ + . . . = log A; in the numerator we have log 2 since the Farey tree is 
a left-right partition, like the ternary of Cantor. 



